function [init,Upnew,Lownew] = Initialization1(pop,dim,Up,Low)
    [M1,N1] = size(Up);
    [M2,N2] = size(Low);
    if N1 == 1
        Upnew = Up*ones(1,dim);
    else
        Upnew = Up;
    end
    if N2 == 1
        Lownew = Low*ones(1,dim);
        else
        Lownew = Low;
    end
    [zonghe] = [Lownew;Upnew];
    for i = 1:dim
        a = rand(pop,1)*(zonghe(2,i)-zonghe(1,i))+zonghe(1,i);
        init(:,i) = [a];
    end
end
